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ABSTRACT 

We have developed a model for the variability of gamma ray emission in jets of active 
galactic nuclei in which the variability arises as a result of photon-photon pair produc- 
tion interactions with X-rays emitted by a hot spot in the inner part of the accretion 
disk. As the hot spot orbits around the central engine, the amount of absorption varies 
periodically. Our model may account for the observed variability of TeV emission from 
Markarian 421 and other blazars detected by the EGRET instrument on the Compton 
Gamma Ray Observatory, as well as correlated variations of X-ray and TeV gamma 
flux from Markarian 421. Quasi-periodic variations observed in Markarian 421 enable 
us to place an upper limit on the black hole mass of (2 — 30) x 10^ M©. 
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1 INTRODUCTION 

Two BL Lac objects, Markarian 421 and Markarian 501, 
have recently been detected at TeV energies (Punch et 
al. 1992, Petry et al. 1996, Quinn et al. 1996). The obser- 
vations of Markarian 421 show that the emission is vari- 
able on different time scales from weeks to days (Kerrick et 
al. 1995, Macomb et al. 1996, Schubnell et al. 1996, Buckley 
et al. 1996), and even variability over a fraction of an hour 
has been observed (Gaidos et al. 1996). In multiwavelength 
observations of Markarian 421, simultaneous variations in 
X-rays and TeV 7-rays have been reported (Takahashi et 
al. 1996a, b, Buckley et al. 1996). Also, the X-ray flux shows 
quasi-periodic variability on a time scale of ~1 day super- 
imposed on a ~ 1 week time scale decline of the emission 
(Takahashi et al. 1996a). 

This last feature is not easily explained in the mod- 
els which assume that X-rays and TeV 7-rays come from 
a relativistic blob moving along the jet. However, it could 
be explained by a model in which the quasi-periodicity is 
caused by the rotation of a massive star which enters the jet 
and creates a shock at some distance from its surface, as re- 
cently proposed by Bednarek & Protheroe (1996). Another 
natural explanation is possible if the X-rays originate in a 
relatively small hot spot rotating on the surface of the in- 
ner accretion disk. If the inner disk thickness increases with 
radius, quasi-periodic modulation of X-ray emission may be 
observed as the hot spot orbits around the black hole with 
the disk, and is viewed at different inclination angles. The 
hot spot emission may be also intrinsicly coUimated by the 
hot spot's geometry (e.g. the walls of the hot cavity in the 



disk), or by an ordered magnetic fleld, and this may increase 
the level of modulation. 

Since the TeV emission can not originate in the hot spot 
itself (because of strong absorption), it must come from the 
jet above the accretion disk. However the processes respon- 
sible for the X-ray and TeV 7-ray variable emission can still 
be physically linked. It is very likely that significant changes 
in the jet conditions are initiated by changes in the accretion 
conditions at the inner disk. For example, a sudden release of 
energy in the inner disk region can change drastically the vis- 
cosity parameter, propelling huge amounts of accreting mat- 
ter, and frozen-in magnetic field, into the black hole and/or 
the jet. Such an energy release may be caused, for example 
by an explosive reconnection of the magnetic field in the in- 
ner part of the disk (e.g. Haswell, Tajima & Sakai 1992), or 
as a result of a thermonuclear explosion on the surface of a 
compact object (white dwarf, neutron star) which has been 
captured by the accretion disk. If this happens, the compact 
object can accumulate enough matter on its surface by ac- 
cretion in the dense medium to initialize a thermonuclear 
outburst similar to those in scenarios proposed for nova or 
Type I supernovae explosions (e.g. MacDonald 1983, Kovetz 
& Prialnik 1994) . The phenomena mentioned above can cre- 
ate a hot spot (or a hot hole) on the surface of the disk. The 
energetic particles in such a hot spot, bathed in the magnetic 
and radiation field, may emit X-ray photons with the spec- 
trum observed during an outburst in Markarian 421. The 
possible occurrence of low photon energy outbursts in the 
inner accretion disk has been recently proposed by Bottcher 
& Dermer (1995) as a method for localizing the emission 
region of VHE 7-rays in active galactic nuclei (AGN) or 
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Figure 1. (a) Simple picture of X-ray emitting hot spot and emission region of VHE 7-rays located at a distance x from the hot spot. 
The 7-rays, injected at an angle ip to the direction from the hot spot, interact with the X-ray radiation from the hot spot, (b) Schematic 
representation (not to scale) of the central region of AGN. The X-ray emitting hot spot orbits with the accretion disk at a distance Ts 
from a massive blaclc hole. The 7-rays emitted at angle a to the jet axis and at different distances / along the jet can interact with X-rays 
from the hot spot. The location of the hot spot on the disk is defined by its phase </>. 



obtaining an upper limit on the column density of matter 
around the disk if the low energy photons are scattered by 
this matter. 

The TeV emission could be produced in the jet by 

any of the mechanisms suggested as being responsible for 
7-ray production in blazars (e.g. Maraschi, Ghisellini & 
Celotti 1992, Mannheim & Biermann 1992, Dermer & 
Schlickeiser 1993, Sikora, Begelman & Rees 1994, Bland- 
ford & Levinson 1995, Bednarek, Kirk & Mastichiadis 1996, 
Bloom & Marscher 1996, Protheroe 1997). If the TeV 7-rays 
are not produced too far away from the hot spot, then they 
can be absorbed by photon-photon pair production with X- 
ray photons from the hot spot, with the escape probability 
depending on emission angle and energy. The absorption 
strongly depends on the location of the TeV emission region 
relative to the hot spot on the disk and the observer. 

In this paper we discuss the effects of such possible se- 
lective absorption on the escape of very high energy (VHE) 
7-rays, assuming a simple hot spot jet geometry. We start in 
Sect. 2 with the simplest possible scenario: a relatively small 
X-ray emitting hot spot and an emission region at some dis- 
tance above the hot spot which injects TeV 7-rays at some 
angle. We then discuss in Sect. 3 a more complicated, but 
probably more realistic, scenario in which the X-ray hot spot 
is created somewhere in the inner disk and orbits with the 
disk around the central black hole. Here, the TeV emission 
region extends along the jet axis, and injects TeV 7-rays at 
some angle to the jet axis. The 7 ray flux seen by a distant 
observer will then be modulated with the orbital period of 
the hot spot by selective absorption, unless the absorption 
of TeV 7-rays in the accretion disk radiation dominates. We 
compute the optical depth for TeV 7-rays in the accretion 
disk radiation, and put the limits on the inner disk tempera^ 
ture, in Sect. 4. The variability of TeV emission expected in 



such scenarios, and its possible relevance to 7-ray modular- 

tion in compact sources, e.g., VHE 7-ray blazars or galactic 
superluminal sources, is discussed in Sect. 5. 



2 OPTICAL DEPTH OF GAMMA-RAYS IN 
THE X-RAYS EMITTED BY THE HOT SPOT 

Let us assume that the X-rays are produced in a small hot 
spot, and the TeV 7-rays are injected from the blob at a 
distance x above the hot spot, and propagate at an angle 
tp measured from the direction defined by the hot spot and 
the blob (see Fig. la). To obtain some idea about the pos- 
sible level of absorption of TeV 7-rays in a real source, we 
assume that the hot spot emits X-rays with a broken power 
law spectrum (phot cm^^ s^^ keV~^), with an exponential 
cut-off below ^ 20 keV, similar to that observed during an 
outburst in Markarian 421 (Takahashi et al. 1996a,b, Ma- 
comb et al. 1996, Buckley et al. 1996): 

p . ^ f 0.187e-l■^ e<1.65keV; , . 

\ 0.229e-2-3 exp (e/20 keV), e > 1.65 keV. 

The optical depth for a 7-ray photon with energy E~f , mea- 
sured from the point of injection to infinity, as a function of 
the injection angle ^, is given by 

r{E^,iP) = ^J (1-^»^X) ^^ J F^(^,)^^^(^,)de, (2) 

where da f» 180 Mpc is the distance to Markarian 421, D = 
+ + 2H cos^l))-^/^ , cosx ^ {H + cos'i(j)/D, H is the 
7-ray propagation distance in urnts of x, s — 2ei5-y (1 — cos x) 
is the center of momentum frame energy squared, Fx{t) is 
given by Eq. (1), and C7-y(s) is the cross section for photon- 
photon pair production (Jauch & Rorlich 1980). 
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The results for three different 7-ray energies are shown 
in Fig. 2. It is clear that, for reasonable hot spot-blob dis- 
tances, the interaction of TeV 7-rays with hot spot X-rays 
can be important. For example, if a; = lO'^* cm 7-ray pho- 
tons with energy = 1 TeV can escape only within a cone 
with opening angle ~ 2°. For x = 10^^ cm, this escape 
cone has opening angle ~ 17°. Hence, if the X-rays are pro- 
duced in a relatively small region on the accretion disk, then 
the escape of TeV 7-rays may be strongly dependent on the 
geometry of the system. In the next section we discuss a pos- 
sible geometry for the emission regions of X-rays and VHE 
7-rays which may occur in the central regions of AGN. 



3 HOT SPOT ORBITING AROUND THE 
CENTRAL ENGINE 

We now consider a more complicated scenario. Let us assume 
that there is a small region (hot spot) in the inner disk at a 
distance from the central black hole orbiting with angular 
velocity lOs (see Fig. lb) . We neglect the hot spot dimension 
which is assumed to be much smaller than the characteristic 
distance scale, Vs- The angular velocity of the hot spot is 
determined by Vs and the black hole mass, and is given by 
Kepler's law 



= {GMf'\:'''^ = c(^,/2)l/v^^''^ 



(3) 



where c is the velocity of light, and rg = 2GM/c? is the 
Schwarzschild radius of a black hole with a mass M. 

As mentioned above, we assume that the 7-rays are con- 
tinuously injected from the stream of plasma moving along 
the jet, and a is the viewing angle with respect to the jet 
axis (Fig. lb). They may also originate in a blob, or a se- 
quence of blobs, moving along the jet. The TeV 7-rays move 
through the X-ray radiation field of the hot spot located on 
the disk surface at phase angle (j) which depends on time. 



<j)= 4)0+ UJst, 



(4) 



where (po is the initial phase of the hot spot (see Fig. lb), 
and we neglect any inward radial motion of the hot spot. 
As discussed below, this leads to a periodic variation of the 
angle between the 7-rays propagating towards the observer 
and the X-rays from the hot spot leading to periodic vari- 
ation of the 7-ray flux. Kepler's law (Eq. ^ gives us the 
radius of the hot spot in terms of the black hole mass and 
the angular velocity. 



which gives 



rs = 1.4 X lO^A^y^ cm, 



(5) 



(6) 



for an orbital period of ~ 1 day, i.e. corresponding to the 
observed quasi-periodic variability. However, this hot spot 
orbital radius must be at least that of the minimum inner 
accretion disk radius, ~ (0.5 — 3)rg, for maximally rotating 
Kerr and Schwarzschild black holes respectively, giving 



> (1.5 - 9) X 10"M8 



(7) 



Hence we obtain a strict upper limits to the mass of the black 
hole in Markarian 421 of M < 2 x 10* M© (Schwarzschild 
black hole) or M < 3 x 10^ Mq (Kerr black hole). 
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Figure 2. The optical depth for the VHE 7-ray, multiplied by the 
distance x of photon injection from the hot spot, versus the angle 
of photon injection, ip. Separate curves correspond to different 
7-ray energies: 0.32 TeV (dotted line), 1 TeV (full line), and 3.2 
TeV (dashed line). 
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Figure 3. The optical depth, multiplied by the radius of the hot 
spot orbit rs , for 1 TeV 7-rays as a function of distance / of their 
injection point from the black hole, for angles of photon injection 
a = 2°, 5°, 8°, and 13° degrees. Phase (f> = 0° is shown by the 
full lines and at </> = 180° by the dashed lines. 



The optical depth for 7-ray photons can be calculated 

o: 

ds /" (1 - cosx) 



D2 



dH / Fx{e)a^^{s)de, (8) 



where it is convenient to express all distances in units of 
H = h/rg, L — l/rg. The values in Eq. (H) are given by. 



D = {1 + H'^ + L'^ + 2HLsmacos(l))-^^'^, 
and 

cos X ~ {H + L cos a — sin a cos 0) /D. 



(9) 
(10) 
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Figure 4. Transmission probability as a function of pliasc </< of the hot spot. The 7-rays are iiijeeted at a distance I and at angle a. 

(a) Dependence of transmission probability on the hot spot phase for 1 TeV 7-rays, projjagating at a = 7°, for Ts = 10^"^ cm, and for 
log(//rs) = 1 (dotted line), 

TeV (dotted line), 1 TeV (full line), and 3.2 TeV (dashed line). 



vith 



10^ 



, Q = 7°, and log(i/rs) = 1.5. (c) Dependence on 
the injection ang le: a = 5° (dotted line), 7° (full line), and 9° (dashed line), with = I TeV, = lO" cm, and \og{l/ra) = 1.5. (d) 
Dependence on the orbital radius of the hot spot: rs = 3 X 10^^ cm (dotted line), 10^* cm (full line), 3 X 10^* cm (dashed line), with 
= 1 TeV, a = 7°, and log(Z/rs) = 1.5 



In order to gain an impression about the possible level 
of modulation of the VHE 7-ray (^mission, we compute the 
optical depth &t = 1 TeV for two critical locations of the 
hot spot in its orbit around the black hole: </> = 0° , and 180° . 
The results are shown in Fig. 3 as a function of the distance 
of the 7-ray injection point from the black hole I, for four 
different angles of observation a. If the hot spot is located at 
phase = 0°, and a is large, then the optical depth drops off 
monotonically with distance. The dependence of the optical 
depth on /, for when the hot spot is located at (j) = 180°, 
shows characteristic cusps at distances which correspond to 
the 7-ray photons and the X-ray photons from the hot spot 
moving such that the angle between their directions is very 
small. The location of the cusps shifts to smaller I with in- 
creasing a. The optical depths for = 0° and = 180° 
converge at large distances I from the black hole because 



the angles between interacting photons become comparable. 
The largest modulation of the 7-ray signal is expected at dis- 
tances for which the differences between the full and dashed 
lines, computed for the same value of a, are highest, and 
have r > 1 for = 0° (full Unes) and r < 1 for = 180° 
(dashed lines). 

Closer inspection of Fig. 3 allows us to choose the pa- 
rameters in this scenario (rs, I, a) for which absorption ef- 
fects can significantly modulate the 1 TeV 7-ray signal with 
the orbital period of the hot spot. In order to show the 
level of possible modulation, we compute the transmission 
probability for 7-ray photons (e~^), selecting as a basic ex- 
ample the case where E-, = 1 TeV, log(//rs) = 1.5, and 
rs = lO^** cm. We investigate the dependence of the trans- 
mission probability as a function of phase (f> for different 
values of the parameters mentioned above. Fig. 4a shows 
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the level of modulation for 1 TcV 7-rays injected at angle 
Q = 7°, for three distances from the black hole: log(?/rs) = 
1.0, 1.5, and 2.0. At / = lOr^ (dotted hne in Fig. 4a), the 
modulation is very strong and, in fact, the 1 TeV 7-rays can 
only escape when the hot spot is close to phase <j} — 180°. 
At larger distances the level of modulation is lower (a fac- 
tor of ~ 15 for I « 30rs; full line in Fig. 4a), and becomes 
negligible at I = lOOrs (dashed line). 

The dependence of the transmission probability on the 
7-ray energy is shown in Fig. 4b. The transmission probabil- 
ity for low energy 7-rays is high, but the level of modulation 
of the 7-ray signal by the orbital period of the hot spot is 
low. For example, for E., = 0.32 TcV, the modulation is 
only by a factor of ~ 3, in comparison to a factor of ~ 15 
at 1 TeV and much higher at 3.2 TeV. The transmission 
probability depends strongly on the angle of 7-ray injection 
(see curves for a = 5°, 7°, and 9° in Fig. 4c). However, the 
general pattern of modulation of the 7-ray signal with or- 
bital phase is conserved. In Fig. 4d we show the dependence 
of the transmission probability on the orbital radius of the 
hot spot for = 1 TeV, a = 7°, and log(//rs) = 1.5. As 
expected, when the transmission probability is lowest (for a 
more compact source; see dotted line for rs = 3 x 10^'' cm) 
the level of modulation of the 7-ray signal with the orbital 
period is highest. 

One must keep in mind that, in the present scenario, 
the variable signal will be multiplied by the 7-ray emission 
efficiency which is likely to decline with the distance along 
the jet ajcis. However, the effects described above may be 
irrelevant if the radiation from the accretion disk in which 
the hot spot is located, is too intense. In the next section we 
shall discuss possible absorption of the TeV 7-ray signal in 
the accretion disk radiation. 



4 ABSORPTION OF VHE GAMMA RAYS IN 
THE DISK RADIATION 

Within the central regions of the AGN, interactions with 
photons directly from the accretion disk can be important 
(Becker and Kafatos 1995, Bednarek 1996; sec Bcdnarck 
1993 for an earlier discussion of 7-ray escape from the radi- 
ation field of an accretion disk surrounding a neutron star 
in an X-ray binary source). In order to check for which ac- 
cretion disk parameters the modulation effects of VHE 7- 
ray emission by the X-rays from the hot spot, discussed in 
Sect. 3, are not swamped by absorption in the accretion disk 
radiation, we compute the optical depth, t, for 7-rays in 
the disk radiation (e.g. Bednarek 1993, Protheroe & Bier- 
mann 1996). We assume the surface of the accretion disk 
emits black body radiation with temperature varying with 
radius as T(r) = Ti„{r/ri„)~^^*, where Tin is the tempera- 
ture at the inner disk radius Vin (Shakura & Sunyaev 1973). 
The VHE 7-rays can escape if they are injected at distances 
farther than the so-called the "radius of the 7-sphere", r^, 
which is defined by the following condition (e.g. Blandford 
& Levinson 1995, Bednarek 1996), 

r{E^,a,r^,)^l, (11) 

where t is the optical depth to photon-photon pair produc- 
tion for 7-rays propagating at angle a, measured from the 
disk axis, through the radiation field of the accretion disk. 



The radius of the 7-ray photosphere, expressed in units 
of Tin, is shown in Fig. 5 as a function of 7-ray energy for 
different inner disk temperatures and radii. The angle of 7- 
ray injection has been fixed at a = 7° because our results for 
the modulation of the 7-ray signal presented in Figs. 4(a)- 
(d) have been shown mainly for this angle. As expected, 
the radius of 7-ray photosphere increases with 7-ray energy, 
and also with intensity of the radiation field of the disk (i.e. 
higher Tin and/or ri„). For the modulation of VHE 7-rays 
by interaction with X-rays from an accretion disk hot spot 
to take place the emission point of VHE 7-rays must be 
outside the 7-ray photosphere, i.e I > r~i{E.y). 



5 CONCLUSION 

We have shown that if the X-ray emission observed from 
Markarian 421 during the outburst stage (Takahashi et 
al. 1996a,b) originates on the surface of the accretion disk, 
then it creates a sufficiently strong radiation field for the 
absorption of VHE 7-rays provided they originate in the jet 
near the accretion disk (Sect. 2). We have discussed a simple 
geometrical model in which a small hot spot emitting X-rays 
orbits with the accretion disk around the central black hole, 
and the VHE 7-rays are injected at some distance along the 
jet at a small angle to its axis. We have shown that absorp- 
tion effects on VHE 7-rays in the radiation field of a hot 
spot, with an X-ray spectrum similar to that of Markarian 
421, can become very important for some sets of model pa- 
rameters, i.e. angle of injection a, location of injection point 
/, and radius of the orbit of the hot spot (see Fig. 3). Ab- 
sorption becomes stronger for higher energy 7-rays and may 
result in a change in the shape of the VHE 7-ray spectrum. 

Modulation of the 7-ray emission with the orbital pe- 
riod of the hot spot is also very significant, and may result 
in simultaneous quasi-periodic oscillations of X-ray and TeV 
7-ray emission. Note that ~ 1 day quasi-periodic variabil- 
ity has been detected by ASCA in X-rays from Markarian 
421 (Takahashi et al. 1996a), and enables us to place an up- 
per limit on the black hole mass for the present scenario of 
M <2x 10** M© (Schwarzschild black hole) ot M <3x 10^ 
M0 (Kerr black hole). Modulation effects such as those de- 
scribed in the present work are expected provided the disk 
radiation itself does not produce a sufficiently intense radi- 
ation field for VHE 7-ray absorption. For the parameters 
discussed above, the disk radiation is not important if the 
disk inner temperature is less than ~ 2 x 10** K for a disk in- 
ner radius rin = W^'^ cm. The direct accretion disk radiation 
is observed in some quasars e.g., in quasar 3C 273, which is 
a 7-ray emitting active galaxy and has an inner disk temper- 
ature (2-3) X lO^K (Shields 1978, Malkan & Sargent 1982). 
The disk radiation fields of the BL Lacertae objects, e.g. 
Markarian 421, are significantly weaker since their spectra 
arc continuous, without strong evidence of emission lines or 
other thermal features. 

The TeV 7-rays may also be absorbed on radiation scat- 
tered by matter distributed around the accretion disk (e.g. 
Sikora, Bcgclman & Rees 1994, Botcher & Dermer 1995). 
However, these radiation fields, if important, arc difficult to 
consider in detail because they depend on parameters which 
are not known precisely. In some sources, where the bulk of 
the infrared emission comes from a dusty molecular torus 
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Figure 5. Radius of the 7-ray photosphere r-y, expressed in 
units of the disk inner radius ri„, as a function of the energy of 
7-ray photons. The disk is a Shakura-Sunyaev type (Shakura & 
Sunyaev 1973) with the temperature Tj„ at the inner ra<lius Tin ■ 
Separate curves correspond to different temperatures at the disk 
inner radius: Ti„ = 2 X 10* K, and ri„ = 10^* cm (full curve); 
Tj„ = 2 X 10* K, and ri„ = lO^^ cm (dotted curve); Tj„ = 10* 
K, and ri„ = 10^* cm (dashed curve). 



surrounding the central engine, VHE emission from points 
below the top of the torus would be absorbed by photon- 
photon pair production on infrared radiation from the torus 
(Protheroe and Bicrmann 1996). However, in the absence 
of such a torus, or if the temperature of the torus is low, 
modulation of VHE 7-rays by interaction with X-rays from 
an accretion disk hot spot as described in the present paper 
should be possible. 

The stream of relativistic jet plasma, or sequence of 
blobs, from which the VHE 7-rays might originate does not 
need to move rectilinearly along the jet axis. It may emerge 
from the region of the hot spot and move outwards from the 
disk following, e.g., helical magnetic field lines. Farther from 
the disk, the angular rotation of the plasma stream can be 
slower than that of the inner disk. In such a situation the 
period of the hot spot's orbit, responsible for the modulation 
of observed X-ray flux in this scenario, may be shorter than 
the period of the TeV 7-ray modulation. 

An interesting feature of the present scenario is that it 
predicts a time delay of the peak TeV 7-ray flux with respect 
to the peak X-ray flux. This is caused by the different paths 
to the observer from the X-ray and 7-ray emission regions. 
The delay depends on the distance of the 7-ray emitting 
region from the disk I, the radius of the hot spot's orbit Vs, 
and the angle of observation a, and is given by 
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At = {r'i + Z^)'/^[l - cos(6l - a)]/c, 



(12) 



where, tanO — Ts/I. For example, if = 10^* cm, a — 7° 
and I = 30rs , then the expected delay of the maximum in 
7-ray flux is At « 40 minutes, and for I = IQvs the delay is 
At w 6 minutes. 
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